hep-ph/9t)lU44U| (revised March 7 1997) 
Breaking Flavour Symmetry Spontaneously. 



Nils A. Tornqvist 

Research Institute for High Energy Physics, SEFT FOB 9, FIN-OOOI4, University of Helsinki, 

Finland 

: (7 March 1997) 



Abstract 

A new mechanism for spontaneous breaking of flavour symmetry is demon- 
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. strated. An exactly flavour symmetric model with degenerate bare nonets 

. and with sufficiently strong tri-linear meson couplings is shown to lead to self- 

. consistency equations which are unstable. Instead there exists a stable solu- 



tion, which break flavour symmetry spontaneously in the mass spectrum. For 
Oh. a C-degenerate meson spectrum the stable mass spectrum obeys the Okubo- 

Zweig-Iizuka (OZI) rule and the approximate equal spacing rule. 
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1. Introduction and the NJL model. Any solution of the hadron mass spectrum, con- 
sistent with QCD, must satisfy the constraints of: (z) unitarity, {ii) analyticity, and in the 
massless quark limit also {iii) chiral symmetry and {iv) flavour symmetry. 

Conventionally one breaks {iii-iv) by adding effective non-degenerate quark masses to 
the Lagrangian, wherby the pseudoscalars obtain (small) masses and the degeneracy of 
all flavour multiplets is split. The gluon anomaly contributes to the flavourless sector in 
breaking the symmetry to approximate SU3f, giving especially the rj mesons part of their 
masses. Most of the chiral quark masses are assumed to come from a short distance regime, 
where weak interactions are relevant. Phenomeno logically fermion masses can be generated 
by suitable ad hoc couplings to a Higgs field. Here I shall discuss an alternative way of 
breaking the flavour symmetry. 
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From deep inelastic scattering, including the spin problem of the proton, we know that 
mesons and baryons are not simple qq and qqq quark model bound states, but have, in 
addition to the valence quarks, large components in the form of a multiquark qq sea and 
gluons. Thus a constituent quark of the naive quark model must, in fact, be composed of 
many quarks, antiquarks and gluons. This requires a multichannel formalism, which couples 
qq and multiquark states. 

Recently it was shown that one can obtain a good understanding of the light scalar meson 
nonet when one uses the constraints {i-iv) above provided one includes in a coupled chan- 
nel framework all light two pseudoscalar thresholds with flavour symmetry broken mainly by 
the pseudoscalar masses. Similar models for other multiplets also improve our understanding 
of the hadron spectrum, even for heavy quarks. One good example is the anomalously large 
T(5)-T(4) mass splitting 0], where single channel potential models clearly fail, because of 
the opening of the BB, BB* etc. thresholds. 

In such models one needs a smaller bare — quark mass term than usual, and 
this becomes smaller the more hadronic thresholds are included. With sufficient number 
of higher thresholds (obeying the empirically approximate equal spacing rule) the SU3f 
breaking usually attributed to nis — rrid could be generated by hadronic loops. This opens 
up the possibility that the bare rris — rriii could be put equal zero, and that flavour symme- 
try breaking could be generated spontaneously, through the self-consistency equations for 
hadron propagators, which include hadronic thresholds. 

The well known Nambu-Jona-Lasinio (NJL) model is a standard reference model for 
for chiral symmetry breaking and for dynamical symmetry breaking in general. Although 
I shall not use this model here the instability mechanism in the original NJL model has 
similarities with the mechanism I shall discuss, and therefore I recall here some well known 
results from the NJL model. For reviews see Ref. 0. With a small chiral mass (/x) the well 
known gap equation for the quark constituent mass (M) is 

Here A is the cutoff and onjl measures the strength of the four fermion interaction. The 
function F takes into account the dependence on the fermion mass in the tadpole-like loop. 
It will be defined below in a dispersive framework. 

Eq. ([1|) has for unjl > 1, = 0, apart from the trivial, unstable solution M = 0, a 
massive and stable solution. This is most easily seen by plotting the left side of the equation. 
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MoutI against M(„/A of the right side as in Fig. 1. For aNjL > 1 the slope of the curve at 
the origin is > 1. Then starting at any value for Mj„/A and iterating the equation one ends 
up at the massive solution as shown in Fig. 1. On the other hand if a jvjl < 1 one obtains 
the dashed curve and the trivial massless solution is the stable one. 

At the same time as the fermions obtain mass, poles are generated in the pseudoscalar 
and scalar qq channels. The pion is massless as required by chiral symmetry while the a 
has a mass of 2M. In generalizations one expects a massless O""*" multiplet and massive but 
degenerate multiplets of other J^'-' when the SU{Nf) remains unbroken. In the following I 
shall not rely further on the NJL model. 

2. A flavour symmetric scalar model. We now construct a meson model, where all 
flavour related bare qiqj states are degenerate, and which couple to each others through 
flavour symmetric couplings. These mesons can then be described in the ideally mixed 
reference frame, i.e. the flavourless states are simply unmixed (like ss, cc). Consider 
now meson loops as shown in Fig. 2 and assume that disconnected quark-line loop diagrams 
which violate the OZI rule can be neglected. To be exact this requires that we have at 
least two degenerate nonets with opposite charge conjugation, and equal F-like and D-like 
flavour couplings as described in more detail in Ref Otherwise singlet states are shifted 
differently from nonsinglet states and quark-line disconnected loops (like (p —>■ u through 
strange intermediate states) do not vanish. Here this C-degeneracy is assumed for simplicity, 
but it is relaxed in Ref. [Q for a model with three flavours, where the OZI rule is broken, 
but the isospin subgroup remains unbroken. 

Taking into account these loops, or the vacuum polarization diagrams, the meson spec- 
trum should still be consistent with being degenerate, if one disregards possibility of insta- 
bility. This is satisfied, since also the thresholds rriik + are degenerate in flavour, and 
when summed over k, each meson ij (using the shorthand qiqj = ij) gets an equal contri- 
bution from the vacuum polarization diagrams. Of course to have this result, not only the 
couplings but also any cutoff. A, or subtraction constant involved must be independent of 
flavour. Then the expression for the self energy, which include quantum loops can be the 
same for all ij, and also the renormalized masses (which now include "unitarity effects") can 
remain degenerate. Consequently there is one very symmetric situation, where one knows 
the solution, and our equations must allow for this trivial self-consistent solution. But is 
this solution stable? 1 shall show that for sufficiently large coupling it is not! 

3. A new mechanism for spontaneous symmetry breaking. The simplest, model is ob- 
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tained for the case of only two flavours and scalar mesons coupling to two-body thresh- 
olds of scalar mesons. Denote the two flavours by 1 and 2. Then there are four mesons, 
ij = 11, 12, 21 and 22, and for each meson there are meson-meson thresholds with different 
quark content, rriik + mkj, k = 1,2. The inverse propagators, get contributions from 

the meson loops due to these thresholds. One finds a sum of two contributions (cf. Fig. 2) 

P^:^\s)=ml~s + ^ ^ F{s,ml,ml.,A) , (2) 

fe=i 

where g is the coupling for each threshold, and mo is the common bare mass. The function 
F must have the unitarity cut, which is proportional to two-body phase space. A simple 
model for F, which I use for the demonstration in this paper, has an imaginary part given by 
s X k / y/sN [s] , where k/ y/s is the usual two-body phase space factor, and N{s) (in the N/D 
formalism) is for simplicity, N{s) = — /c as in an effevtive theory, and a real part given by 
a dispersion relation. Denoting as usual the Kallen function {x — y — z)^ — \.yz by A(x, z), 
the imaginary part is thus: Im[F] = Xi{s,ml,ml)9{\)9{A'^ — A/(4s)). It vanishes below 
threshold (the first 6 function) and for momenta above the cutoff (the second 6 function). 
The real part is determined by the dispersion relation: 

Re F(s,mi,m2, A) = -P/ . (3) 

n J0<A<4s'A2 S' — S 

The function F can be evaluated analytically and is simply related (up to subtractions and 
factors of s in the imaginary part) to the Chew-Mandelstam, and H{s) functions, which also 
appear in the literature. With this F the coupling constant g is dimensionless. 

The physical meson masses (poles) are given by the zeroes of the inverse propagators of 
Eq.(0). Now eliminate the universal bare mass mo by fixing the mass of one of the mesons, 
say mil. One then obtains the following self-consistency equations (cf. Fig. 2): 

g2 Nf=2 

Q = Pij\s = ml)=ml,-s + f- [F{s,ml,ml,A) - F{ml„ml„ml,,A)] . (4) 

k=i 

By construction mn is always the solution of Pi7^(mJ^)= 0. The two other masses, mi2 
and m225 are then determined by the remaining two equations. From the above discussion 
the degenerate solution with equal rriij = niu is one solution. But this need not be the 
only solution, and furthermore the symmetric solution need not be stable. Denoting a small 
variation from the symmetric solution in the threshold masses by 5ij = mfj — rriQ, which 
results in a shift 5°*^* in the pole positions at the right hand side of Eq. (|), self-consistency 
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of course requires 5°"* = 5ij and stability 5°"* < 5ij. I.e., if the latter inequality is satisfied 
then starting from some 5ij 7^ and iterating one converges towards the symmetric solution 
5°"* = 5ij = 0. On the other hand, if 6°j^^ > 6ij the symmetric solution is unstable. Denoting 
by Fg = dF/ds\s=jn2^jn2 and by F^a = dF{s,ml,ml)/dml\g^^2^^2 the self-consistency 
equations can be written 

^Ti^f^^^ - 1] + Et^^fc + ^.fc] = • (5) 

This self-consistency requires that small deviations from the symmetric solution must satisfy 
the equal spacing rule 



C'-C* = ^('^*-C*), (6) 



while the symmetric solution is unstable if 

NfF 



> 1 . (7) 



8u -NfFs + A-K/g^ 

This quantity is plotted in Fig. 3 for large values of the coupling constant. Here r is the 
slope at origin of the function plotted in Fig. 4. Equivalently one can write this condition 
as a bound on (yf^/(47r): 

N^l->[F, + F^2]-K (8) 

The left hand side is always positive for any reasonable function F, with correct threshold be- 
haviour. Thus the instability occurs for sufficiently large coupling g. In it was found that 
typical coupling constants such as gpT^.,^ and (^o-TTTr very well satisfy this instability condition, 
indicating that the instability we discuss actually occurs in Nature. Since this instability 
depends only on the threshold behaviour of relativistic phase space, it is a fundamental prop- 
erty related only to the underlying Lorenz invariance and flavour symmetry of the model. 
When g satisfies the condition one must look for another stable solution. This can only be 
done numerically since the function F is nonlinear already in the simplest possible models. 
To find the stable solution one must solve the nonlinear equations. This can be done by 
iterating the Eqs. (|^), starting with some value off the symmetric solution (See Fig. 4), in a 
way analoguous to Fig. 1. Using the input masses niij^in for the threshold masses in the loop 
one calculates for which s = mfj one has zeroes in the inverse propagators Eq. (^) . The 
latter are then in the next iteration are used as input masses for the thresholds. Of course, 
a direct way is to solve without iterations the nonlinear equations on a computer. These 
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solutions are shown in Fig. 5a for the functon F defined as above and in the hmit when g is 
very large. 

4-. The symmetry is broken also in the wave functions. How can this unsymmetric 
solution arise from a completely fiavour symmetric theory? All four mesons are apart from 
the "bare qq seed" composed of clouds of multi-quark pairs in the form of meson-meson pairs. 
But these clouds can be different! One can write for the wave functions decompositions: 

Kj>oc ^[V^|zl,lj> +V^|z2,2j>] + |«j> . (9) 

The coefficients z]^ are generally different and can be computed from the relative slopes 
with respect to s of the function F evaluated at the stable solution z]^ oc — ^F(s = 
mfpmff^,m\^, K). The normalized probabilities = zl l{z-l + Z2^ for the stable solu- 
tion of the model discussed are shown in Fig. 5b. Only for the unstable, symmetric, solution 
are the wave functions the same for all ij. The equations are also self-consistency equa- 
tions, i.e., the quantities ij on the right hand side stand for a collective multiquark state 
obtained when this quantity is iterated into these quantities on the right hand sides. Thus 
the true physical Fock states have components with an arbitrary number of virtual quarks 
or mesons in their wave functions. 

5. Concluding remarks. The output stable spectrum is ideally mixed and obeys approxi- 
mately the equal spacing rule — mu ^ 2(m2i — mn), as one should expect from Eq. (^. 
Thus just as in the real world one can define the constituent quark mass as approximately 
Mi ^ mjj/2, whereby rriij ^ Mj + Mj. 

A natural very important question which arises in this connection is: Where are the 
Goldstone degrees of freedom and the Goldstone bosons expected whenever a symmetry is 
spontaneously broken? In Ref. ^ I argue within a scalar QCD model that actually the 
scalar or longitudinal confined gluons are the would-be Goldstone bosons, not scalar mesons 
carrying fiavour as is usually expected. 

In the model presented above there were only S-wave thresholds. But the instability 
increases (r grows) if the threshold involve angular momentum factors k'^^, since this in- 
creases the sensitivity to the thresholds and in particular in Eq. ([^ grows. Above only 
scalar mesons were considered. Adding other multiplets of different J^*^ does not change 
the picture qualitatively, since most multiplets obey approximately the equal spacing rule. 
After summing over all thresholds one has a similar behaviour in P~^. Effectively this is 
approximately equivalent to increasing g and A in the present model. 
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The instability occurs for any number of flavours Nf. Of course, how the symmetry group 
is broken down to a lower symmetry depends on details of the function F as one moves off the 
symmetry point, and how the mixing between the flavourless states evolves. In numerical 
experiments with Nf = 31 find that generally an SU2 subgroup remain unbroken (See 
Ref. ii. 

The assumption of C-degeneracy of multiplets which was made above is certainly not 
exactly true in reality. In the real world this is certainly broken, although on the average 
there can be an almost equal contribution from F- and D-coupled flavour thresholds. There- 
fore the OZI rule is approximately valid for most multiplets. I relax the C-degeneracy in 
Ref. using a model for the ground state mesons, including the pion. Then the OZI rule 
must be violated by loop diagrams, but the isospin subroup remains unbroken. 

Since the present mechanism requires smaller quark masses it may be important in the 
resolving of the strong CP problem, for which one resolution would require a vanishing quark 
mass. Also, there need be no contradiction with standard relations between pseudoscalar 
masses and usual quark masses since the usual quark masses are "measured" through 
these relations. Including the loop effects discussed here, the scale of the quark masses 
decreases while B in m'^{0~~^) =B{mqi +mq2) increases. In lattice gauge theory p one also 
calculates light quark masses, but there one either uses a quenched approximation, or make 
crude approximations for quark loops. In the present work the loops play a crucial role, 
without which the mechanism of spontaneous breaking would not work. Therefore there is 
no contradiction with present lattice calculations. 

Finally it should be pointed out that this work does revive old ideas from the 60'ies and 
70'ies, where one hoped to bootstrap the hadron spectrum through similar self-consistency 
equations 0,|lOl, but to my knowledge no demonstration like the one I have presented here 
was ever presented. 

I thank especially Geoffrey F. Chew and Claus Montonen for useful discussions. 
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FIGURES 
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FIG. 1. Graphical solution of the gap equation. For aNJL > 1 the slope of the soHd curve 
at the origin is > 1. Then the stable quark mass M is massive, otherwise one finds the trivial 
massless solution (dashed curve). Here unjl = 1-30 from which M/A = 0.48, which fits /^r = 93 
MeV. Then A =0.653 MeV if M = Mjv/3. 



FIG. 2. The self-consistency equation diagrammatically. Iterating the equation gives a sum 
of multiloop diagrams. 
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FIG. 3. The ratio r at the the symmetric point {x = y = z = mg/A) is shown when g is very 
large. The fact that r > 1 impies unstabihty of flavour symmetry for all A and spontaneous flavour 
symmetry breaking occurs. 
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FIG. 4. The spontaneuos breaking of flavour symmetry for two flavours, when mn/A = 0.5 
and g very large. Compare this figure with Fig. 1. The slope of the curve at the origin is given by 
r, which is always > 1 for all A (Fig. 3) . The symmetric point of equal quark masses is unstable 
(like the massless point in Fig. 1), while the stable solution has unequal masses. 
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FIG. 5. a) The spontaneous splitting of meson masses rriij/K as a function of the hghtest 
meson mass in units of. Note that the stable solution approximately satisfies the equal spacing 
rule, 777,22 — fni2 ~ ?tii2 — fnw, in accord with physical mass splittings. The unstable solution 
with degenerate masses, and same subtraction constant, is shown by the dashed curve, b) The 
normalized probabilities Z]^-' of Eq. (6). 
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